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ABSTRACT 


The  numerical  procedure  of  calculating  the  second- 
order  steady  forces  is  presented.  The  computer  program 
is  developed  on  the  basis  of  Lin  and  Reed's  theory (1976) . 

Ship  motion  and  diffraction  potentials  are  required  as 
input  data  for  Kochin-f unction  calculation.  In  order  to 
avoid  the  irregular  frequencies  which  are  associated  with 
Frank’s  close-fit  method (1967) ,  a  modification  which 
extends  the  source  distribution  onto  the  calm  waterline 
inside  of  a  body  is  made.  For  the  diffraction  problem, 
instead  of  a  Helmholtz  equation,  a  two-dimensional  Laplace's 
equation  is  used.  Numerical  computation  for  the  head-sea 
case  shows  the  same  trend  of  experimental  data  throughout  the 
frequency  ranges,  but  its  magnitude  is  much  larger  than  that 
of  the  experiment. 

ADMINISTRATIVE  INFORMATION 

This  project  was  authorized  by  the  Naval  Sea  System  Command  with 
funding  under  Element  61153N  Project  SR0230101  and  identified  as  Work 
Unit  Number  1524-680, 


INTRODUCTION 


It  is  a  wall  known  fact  that  when  a  ship  is  navigating  in  a  seaway, 
the  engine  power  must  be  increased  by  a  considerable  amount  in  order  to 
keep  the  same  speed  of  advance  as  in  calm  water.  Especially  for  a  high¬ 
speed  vessel  or  a  ship  on  a  tight  schedule,  it  is  important  to  predict  the 
additional  power  needed  in  advance  not  only  for  navigational  safety  but 
also  for  economical  reasons.  The  steady  force  induced  by  waves  in  the 
beamwise  direction  is  the  lateral  drift  force.  One  may  state  that  added 
resistance  is  important  from  the  standpoint  of  powering  and  design,  and 
drift  force  from  the  standpoint  of  seakeeping  and  control.  Both  the  lateral 
drift  force  and  the  added  resistance  arise  from  the  ship  motion  and 
diffraction  of  waves. 

These  two  forces  are  not  seemingly  related  to  each  other,  but  from  a 
mathematical  point  of  view,  they  can  be  analyzed  by  identical  method. 

In  what  follows  we  shall  call  both  added  resistance  and  lateral  drift  force 
as  the  second-order  steady  forces.  Because  of  the  complexity  of  the 
problem,  there  have  been  only  a  few  efforts  to  study  the  second-order 
steady  forces  analytically  or  experimentally  in  the  past.  The  traditional 
practice  accepted  by  naval  architects  has  been  that  power  increase  in  a 
seaway  is  between  15-30%  of  the  power  required  for  calm-water  resistance. 

Fortunately,  the  advent  of  large  computer  facilities  and  the  rapid 
growth  of  new  computational  technique  for  predicting  ship  motion  make  it 
possible  to  calculate  the  second-order  steady  forces  analytically.  In  the 
last  decade  much  work  has  been  devoted  to  solving  this  problem,  nevertheless 
all  the  studies  have  remained  limited  in  scope  and  achievement.  In  1976 
Lin  and  Reed  presented  a  new  approach  for  evaluating  the  second-order 
steady  forces  in  oblique  waves.  The  forces  are  derived  from  linear 
momentum  consideration.  The  second-order  steady  forces  are  obtained  in 
terms  of  the  Kochin  function  H(u,A)  by  taking  a  time-average  of  the 
periodic  forces  and  invoking  the  method  of  stationary  phase  evaluation  of 
the  potentials  at  a  large  distance  from  the  ship.  The  computation  of 
second-order  steady  forces  will  be  based  on  the  formulae  derived  by  Lin 
and  Reed,  for  their  approach  is  not  only  mathematically  sound  but  also 
much  more  versatile  than  any  other  method. 
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The  validity  of  the  computed  results  is  confirmed  by  '-becking  with 
avai’aK^  model-experimental  and  other  theoretical  values,  and,  also,  the 
effectiveness  of  the  computer  program  as  an  inexpensive  way  of  obtaining 
useful  information  is  ascertained.  The  computer  program  developed  here 
may  be  utilized  for  such  practical  application  as  establishing  data  base 
for  the  design  of  a  ship  for  a  given  route  and  sea  state.  The  work  will 
be  presented  in  the  order  of  a  brief  outline  of  the  mathematical  problem, 
description  of  the  procedure  of  numerical  calculation,  presentation  of 
the  numerical  results,  and  finally  a  summarization  of  the  findings  of  the 
study. 

REVIEW  OF  MATHEMATICAL  FORMULATION 

A  brief  outline  of  the  most  important  aspects  of  the  theory  will  be 
presented  here.  More  detailed  derivations  can  be  found  in  the  report  by 
Lin  =>nd  Reed  (1976). 

The  problem  to  be  considered  here  is  that  of  a  ship  moving  at  constant 
forward  speed  U  with  arbitrary  heading  in  a  plane  of  progressive  waves,  as 
illustrated  in  Figure  1.  Main  assumptions  and  restrictions  in  the  theory 
are  listed  below: 

(1)  The  usualy  ideal-fluid  assumption  is  made,  permitting  the  use 
of  potent ial-flow  theory; 

(2)  The  ship  has  small  displacement  from  the  eqtdl ibr ium  position  and 
both  the  incoming  waves  and  those  created  by  the  ship  are  small; 

(3)  The  ship  is  sufficiently  "slender"  so  that  each  section  can  be 
treated  as  a  two-dimensional  "strip"  with  no  interaction  between  them; 

(A)  The  response  of  the  ship  to  the  incident  wave  is  linear. 

Let  Oxyz  be  a  right-handed  rectangular  coordinate  system  translating 
with  the  mean  position  of  the  ship.  The  origin  is  located  on  the  calm 
water  surface.  The  x-axis  points  forward  and  the  z-axis  points  vertically 
upward . 

Let  'l>  be  the  velocity  potential  given  by: 

<J>(x,v,z,t)  =  -Ux  +  <J>(x,y,z,t)  (2.1) 

* 

References  are  listed  on  page  23 
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and  the  absolute  velocity  of  the  fluid  is  described  by: 


q(x,y,z,t)  =  V<J>  =  -Ue^  +  V<p,  (2.2) 

where  V<j)  =  v  =  ue^  +  +  we^  represents  the  perturbation  fluid  velocity, 

and  e^,  and  the  three  unit  vectors  in  Oxyz  frame.  For  convenience, 
the  perturbation  potential  <J>  is  decomposed  into  a  steady  part  <{>s  and  a 
time-dependent  part  <J>T>  i.e., 

4>(x,y,z,t)  =  (J>s(x,y,z)  +  <|>  (x,y,z,t)  (2.3) 

where  <J>  describes  the  disturbance  due  to  the  steady  forward  motion  of 
the  ship  in  calm  water  and  <J>  can  be  further  decomposed  as: 

=  +  +  $M  =  ♦i  +  (2‘4) 

where  is  the  potential  of  the  incoming  plane  progressive  wave,  <p  the 
diffraction  potential  due  to  the  presence  of  the  ship,  and  (J>  the  forced- 
oscillation  potential.  The  incoming  wave  potential  <(>  is  given  by: 

=  Re{  (x,y,z)exp(iwt)  }  ,  (2.5) 

where 

ij;  =  -  exp{kQz  -  ikQ(xcosB  +  ysinB)  }  (2.6) 

2 

in  which  A  is  the  wave  amplitude,  a  the  wave  frequency,  =  a  /g  the  wave 
number,  8  the  wave  heading  angle  counted  from  the  positive  x-axis  in  the 

counterclockwise  direction,  g  the  gravitational  acceleration  and  to  the 
frequency  of  encounter.  Note  that  a  ship  advancing  through  regular 
sinusoidal  wave  will  not  respond  at  the  frequency  of  incoming  wave,  instead 
at  the  frequency  of  encounter  which  is  defined  as: 


oj  *  a  -  kpUcosB  . 


(2.7) 
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Under  the  assumption  of  linear  response,  we  have 


<p  =  Re{  ip  (x,y,z)exp(iwt)  }.  (2.8) 

D  O 

From  the  linear  momentum  consideration  we  can  derive  the  force,  and 
by  taking  the  time  average  of  the  force,  the  second-order  steady  forces 

•ff  ^ 

<AF>  are  obtained  in  terms  of  and  ip^tp  where  indicates  the 

complex  conjugate.  By  repeated  application  of  the  method  of  stationary 
phase  at  a  large  distance  from  the  ship  and  after  lengthy  algebraic 
manipulation,  the  following  results  are  obtained: 


<AF  >  =  <AFib  >  +  <AFbb  > 


(2.9) 


where 


<AFIB  >  = 


(  CsZ$)lm{  H(B’kn>  }» 


(2  JO) 


<AF. 


BB  > 


8tt 


8u 


tt/2  tt-u0  2 it 

</-/  -f  )du 

-tt/2  tt/2  tt+u0 

Aj(u) 

cosu 
^  sinu 

(1+4tcosu) 

Vo  P 

A~(u) 

cosu 

if  +J  )  du- 

0  TT+Uq 

(1+4tcosu) 

^  sinu 

■  *2'  1 

(2.11) 


in  which 


,  t<  1/4 


u0  = 


-1 


cos  (1/4t)  ,  x>  1/4 


cos  u  and  cos  3  correspond  to  the  added  resistance  and  sin  u  and  sin  3 
correspond  to  the  lateral  drift  force.  The  Kochin  function,  H(u,X), 
and  A ^ ,  A  will  be  explained  in  the  next  section  in  detail. 

As  shown  in  Equations  (2.10)  and  (2.11)  the  forces  are  the  second-order 
quantities,  because  according  to  the  usual  formulation  of  the  ship  motion 
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problem  the  velocity  potentials  are  considered  as  the  first  order  of  the 
wave  amplitude  and  so  is  the  Kochin  function  (sec  Equation  (3.1)). 


KOCH IN  FUNCTION 

The  function  H(u,A)  known  as  the  Kochin  function  is  defined  as 


follows: 


H(U, a)  =jf/'(^°)  -  y*0)  |j-)  exP^pzo+iXp(xOcosu+yoslnu)}ds 


(3.1) 


Here  the  integration  is  over  the  mean  position  of  the  hull  surface  S.,  , 

BO 


Ap  =  - ~ — 2 —  ^  +  2lcosu  “  (1+4tcosu)1/2},  p=l,2 

2x  cos  u 


(3.2) 


where  V  =  to  /g,  r  =  Uw/g,  and  the  plus  and  minus  signs  are  for  p=l  and 
p=2,  respectively.  Figure  2  shows  the  behavior  of  A^  and  A^  as  a  function 
of  u. 

Let  us  first  consider  the  ship-generated  potential  in  detail. 

Using  the  principle  of  linear  superposition  we  can  expand  and  hence 

D 

the  Kochin  function  as  follows: 


<Pn  =  2  >Kv 

B  i=l  1  1  D 


(3.3) 


where  is  the  displacement  of  the  body  due  to  the  motion  in  the  ith  mode, 
and  i=l , 2 , 3 ,4 , 5, 6  represent  surge,  sway,  heave,  roll,  pitch  and  yaw, 
respectively  and 


H(u , A)  =  E  a.H  +  H  . 
i=l  1  1 


(3.4) 


The  potential  iJj  when  multiplied  by  the  complex  amplitude  represents 
the  fluid  disturbance  due  to  the  body  oscillation  in  the  ith  mode,  and  so 
does  the  Kochin  function.  We  can  easily  note  from  Equations  (2. 10) , (2. 11) 
and  (3.1)  that  it  is  very  important  to  predict  the  ship  motion  and  diffrac¬ 
tion  potential  accurately,  because  the  second-order  forces  are  expressed 
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in  terms  of  the  Kochin  function,  and  the  Kochin  function  itself  could  be 
obtained  by  the  integration  of  \ p  and  V  over  the  mean  body  surface  Snrt. 
Therefore,  small  errors  in  the  prediction  of  ship  motion  and  the  potentials 
and  could  result  in  large  errors  in  the  forces. 

SHIP  MOTION  PROBLEM 

Our  ship  motion  computation  has  been  developed  on  the  basis  of  the 
method  of  Salvesen  et  al.(1970).  A  brief  review  of  the  typical 
assumptions  which  are  required  for  the  theoretical  justification  of  the 
strip  theory,  will  be  given  first. 

The  hull  is  assumed  to  be  long  and  slender  such  that  if  B  is  a  typical 
transverse  dimension  and  L  is  the  length,  then 

B/L  =  0(e),  (3.5) 


and 


nx  «  0(e),  n2  =  0(1)  and  n 3  =  0(1)  (3.6) 

where  n^,  n^,  and  n^  are  components  of  unit  normal  vector  on  the  immersed 
ship  surface.  It  follows  that  in  the  neighborhood  of  the  hull 

hr=0<'e)’  |p“  -  0(1)  and  -  0(1)  .  (3.7) 

In  addition  we  shall  make  a  crucial  assumption  that  the  frequency  of 

g 

encounter  is  high,  u  »U  which  means  that  the  wave  length  is  approxi¬ 
mately  of  the  same  order  of  B.  All  these  assumptions  could  only  be 
justified  by  comparing  solutions  with  experimental  data.  We  will  summarize 
the  above  assumptions  in  the  following  manner. 


Original 
Field  equation 


With  Assumption 


dx  9zZ  1 

2  2 

<  H +  )  hm0 

9y  dz  1 

(3.8) 

Free  surface  condition 

(J-UM  |j  -  U^^-g  fj  V0 

dX 

-  s h h  =  0 

(3.9) 

Body  boundary  condition 

3T  =  iwni  +  Umi 

W  “  ^t+^i 

(3.10) 

Radiation  condition 

1/2  *^i 

Lim  R1^  -  ik0^)  -  0 

R-*» 

Lim  (  3-  i  iVi>  *  ° 
y^+» 

(3.11) 

7  t  1/2 

where  R  =  (x  +  y~)  and  n3  is  defined  by 
(ni,n2,n3)  =  n  and  (n4,n5,n6)  =  rxn 


with  n  the  unit  normal  vector  which  is  directed  into  the  body  and  r  the 
position  vector  with  respect  to  the  origin  of  the  coordinate  system  and 
where  m..  *  0  for  i  =  1,2, 3, 4  while  m5  =  n3  and  m6  =  -n2.  is  the 

2-dimensional  generalized  normal  in  the  y-z  plane.  We  note  that  with  the 
above  assumptions,  does  not  depend  on  the  x-coordinate  explicitly. 
Introducing  the  assumptions  into  Equation  (3.1)  and  setting  d^dxdl,  we 
can  also  show  that 
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i 


fix  exp(Uxcosu)/  dl  {  ^  -  (  *,  +  ^  *J  )|jf  } 

L  C(x) 


x  exp{  A(z  +  iysin  u)  },  (3,12) 


in  which  C(x)  is  the  immersed  contour  of  the  cross  section  of  the  ship  at 
station  x,  and  using  body  boundary  condit ion (3 . 10) ,  we  may  obtain 

H^(u,A)  =  f  dx  exp( iAxcosu)^  dl { ( iun^+Um^)  -  (t|^+ 

L  C(x) 

x  A  (N^+iN2sinu) }exp{ A(z+iysinu) ) .  (3.13) 

where  1^=0  for  i=l,2,3,4,  i(;^=  ip^  and  t|^=  -ip^. 

Naturally  in  order  to  calculate  the  Kochin  function  the  velocity 
potential  which  satisfies  Equations  (3.8)  to  (3.11)  must  be  obtained 
beforehand.  There  are  several  methods  in  use  for  finding  solutions  for 
this  boundary-value  problem  such  as  separation  of  variables  and  the  Fourier 
method,  the  method  of  reduction  and  reflection,  the  method  of  Green's 
functions,  the  method  of  multipole  expansion,  and  so  on.  Among  them  the 
method  of  Green's  functions  is  certainly  the  most  flexible  one.  The 
solution  by  this  method  involves  an  integral  equation,  i,e., 

i|/(P)  -  ~J  {i^v(Q)G(P,Q)  -  iMQ)Gv(P,Q)}  ds(Q),  (3.14) 

C(x) 

where  P  *  (y,z),  Q  =  (y0,zQ),  v  the  normal  derivative  in  terms  of  (y0,zQ), 

C(x)  a  contour  bounding  a  two-dimensional  region,  ds  the  arc  length  and 
Green's  function  which  we  denote  G(P,Q)  satisfies  the  following  equations: 

2  2 

(  2  +  2  ^  G(P,Q)  =  6(y-yQ)  6(z-Zq)  ,  (3.15) 

dy  3z 
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(3,16) 


Z  9 

w  G(y,0;y0,z0)  -  g  G  =  0  , 

Lim  {  |^(y’°»yo’Z0)  ±  ik  G  }  =  0.  (3.17) 

dy  u 

y-c+oo 


and  the  solution  is  as  follow  (see  Wehausen  and  Laitone,1960)  : 


G(P, Q)  =  log(r)  -  log  (r') 

f  exP  k0(z+z0)  cos  k0(y_y0)  dk 


-  2PV 


/' 


k-k. 


+  2Tri  *  exp{k0(z+z0)cos(k0(y-y0))} 


(3.18) 


in  which  r  -  {(y-yQ)2  +  (z-Zq)2}'5  ,  r'  «  {(y-yQ)2  +  (z+Zq)2}^,  kQ=u)2/g 

and  PV  indicates  the  Cauchy  principal-value  integral. 

Although  one  can  show  in  many  cases  that  a  solution  to  the  integral 
equation  exists,  a  closed  form  solution  is  usually  not  obtainable,  and 
the  only  possible  way  to  solve  it  is  by  numerical  solution.  With  the 
advent  of  high-speed  computers,  a  numerical  solution  of  integral  equations 
has  become  almost  a  routine  procedure.  However,  one  major  drawback  in 
solving  an  integral  equation  of  the  second  kind  is  the  non-uniqueness 
of  the  solution  when  the  homogeneous  part  of  the  equation  has  nontrivial 
eigen  values.  For  a  body  floating  in  a  free  surface,  John  (1950)  pointed 
out  that  the  integral  equation  involved  admits  non-unique  solutions  at 
the  eigen-frequencies.  He  called  these  eigen-frequencies  irregular 
frequencies.  The  problem  of  irregular  frequencies  has  received  extensive 
investigations.  Paul  Wood  has  demonstrated  by  numerical  computation 
(ref.  Paulling,  1970)  that  the  irregular  frequencies  can  be  removed  by 
extending  the  source  distribution  onto  the  waterline  inside  the  cylinder 
and  imposing  a  rigid  wall  condition  on  it  (see  Figure  3).  Frank  (1967) 
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studied  these  irregular  frequencies  numerically  for  two-dimensional 
cylinders  and  found  that  the  matrix  becomes  ill-conditioned  at  and  within 
a  narrow  band  near  irregular  frequencies.  The  problem  is  important  because 
these  irregular  frequencies  are  not  known  a  priori  for  a  complicated 
geometry.  Ohmatsu  (1975)  has  shown  how  to  avoid  this  difficulty  by 
modifying  the  interior  problem  and  proved  Paul  Wood's  justification  by 
using  Green's  theorem.  Ogilvie  and  Shin  (1978)  have  presented  a  rather 
simple  procedure  that  could  eliminate  irregular  frequencies  by  making 
a  minor  change  in  the  Green's  functions  based  on  a  procedure  suggested 
by  Ur sell  (1953). 

The  reason  we  gave  some  detailed  explanation  for  computation  procedure 
is  that  the  prediction  of  the  second-order  steady  forces  requires  the 
ship  motion  response  and  the  diffraction  potential  as  input  values,  so 
that  the  second-order  steady  forces  predicted  by  a  given  method  may  vary 
considerably  depending  on  the  method  used  for  obtaining  the  motion  and  the 
diffraction  potential  as  also  mentioned  in  Salvesen  (1976).  In  our  compu¬ 
tation  we  adopt  Frank's  close-fit  method  but  we  avoid  irregular  frequencies 
by  adding  a  horizontal  rigid  wall  inside  the  body  by  following  Paul  Wood's 
method.  Figure  4  shows  the  heave  added-mass  coefficients  of  a  circular 
cylinder  with  and  without  a  horizontal  wall.  Figure  5  displays  the  added- 
mass,  damping  coefficients  and  heave,  pitch  magnitudes  for  a  Mariner  hull 
form. 

DIFFRACTION  PROBLEM 

For  the  diffraction  part  of  Kochin  function  Hp(u,A),  we  cannot 
immediately  follow  the  expression  of  Kochin  function  H^(u,A)  for  the 
forced  motion  in  Equation  (3.13).  Before  writing  down  the  desired  final 
form,  let  us  study  the  diffraction  problem  in  some  detail. 

With  the  incident-wave  potential  given  by  Equation  (2.6),  the 
diffraction  potential  is  subject  to  the  condition  that  the  total 
potential  has  zero  normal  velocity  on  the  body  surface.  Since 

the  incident  wave  has  the  factor  exp^ik^xcosB) ,  it  seems  reasonable 
to  expect  that  for  a  slender  ship  and  short  waves  the  diffraction  waves 
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also  have  similarly  oscillatory  behavior  along  the  x-axis.  This  assumption 
is  not  valid  near  the  ends,  but  then  the  assumptions  we  made,  i.e.,  n^ 
changes  slowly  in  the  x-direction  is  not  valid  there  either.  It  is  for 
this  reason  that  we  may  write 


i^D  =  ^jj(y»z) exp(-ikgXcos  6). 


With  this  definition  Y  must  satisfy  the  following  equations: 


V  +  ''Dz,  *  <-‘k0eOS;)'l,)  "  »• 


-  s»Dz  -  o. 


“  "Jn  =  Aw^-iK^sin  P)exp{  k^z-iysin  6)}, 

3Y 

Lim  (  t  ik04u  )  =  0. 

y-*+ou  J 


(3.20) 


(3.21) 

(3.22) 

(3.23) 

(3.24) 


Equation  (J.21)  is  known  as  Helmholtz  equation.  Newman(1970) 
showed  that  the  determination  of  the  sectional  forces  due  to  the  incident 
waves  should  deal  with  a  Helmholtz  equation  in  the  cross  plane  instead 
of  Laplace's  equation  as  the  usual  strip  theory  does.  But  it  is  not 
an  easy  task  to  solve  a  Helmholtz  equation  with  boundary  conditions  (3.22) 
to  (3.24).  Choo  (1975)  solved  the  Helmholtz  equation  and  obtained  the 
diffraction  potential  by  using  an  asymptotic  series  expansion  technique 
for  the  case  of  zero  speed  and  Troesch  (1976)  tried  to  extend  to  that 
forward  motion  at  moderate  speeds  without  obtaining  numerical 
values.  Troesch  compared  his  numerical  computation  with  not  only  experi¬ 
mental  data  but  also  the  solution  of  the  same  boundary  value  problem, 
using  two-dimensional  Laplace's  equation  as  the  governing  equation. 

Figure  6  shows  the  pressure  distribution  for  the  midship  section  of  an 
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ore  carrier  for  L/X  =  1.96  and  6=135°,  45°.  The  integrated  pressure 

forces  are  presented  in  Table  1.  F^  and  Fz  are  the  amplitude  of  the 
sectional  exciting  forces  in  the  horizontal  and  vertical  directions 

respectively  and  are  nondimensionalized  with  respect  to  pgAB/2  where 
B  is  the  sectional  beam. 


Helmholtz 

Laplace 

magnitude 

phase 

(Deg) 

magnitude 

F 

_ 1 _ 

PgAB/2 

1.11 

120 

0.95 

115 

F 

z 

PgAB/2 

0.49 

-55 

0.47 

-57 

TABLE  1  Seccional  Exciting  Force 

In  Figure  7  total  forces  which  are  integrated  over  the  hull  for  the 

o 

Series  60,  CD  =  0.70  are  plotted  for  heading  angle  6=150  and  to(B/2g; 

D 

ranging  from  0.6  to  1.2. 

In  spite  of  the  more  elaborate  numerical  computation  involved  in  the 

solution  of  the  Helmholtz  equation  compared  to  the  solution  of  Laplace's 

equation,  the  results  do  not  seem  to  be  so  different  from  those  of 

Laplace's  equation  as  to  influence  practical  predictions.  Thus,  we  shall 

adopt  Laplace's  equation  in  our  computation,  but  we  have  to  keep  in  mind 

2 

that  neglecting  the(-ikQCOs6)  term  in  Equation  (3.21)  violates  the 
crucial  assumption  we  made  for  the  justification  of  the  strip  theory  in  ship 
motion,  i.e.,  X/L  =  0(e).  Replacement  of  the  Helmholtz  equation  by  the 
Laplace's  equation  saves  considerable  amount  of  computing  time,  because 
once  we  solve  the  forced  motion  problem  numerically,  we  can  immediately 

_ _ _ . _ _ _  i 

*  In  Troesch  8=0(degree)  for  head  seas.  4. 

\ 

.1 
\ 

3 
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obtain  the  diffraction  potentials  by  simply  changing  body  boundary  condition. 
Numerical  procedure  for  this  simplification  will  be  given  in  the  following 
section.  In  passing,  the  Kochin  function  for  the  diffraction  part  will  be 
expressed  as: 


dll>. 


CC  3 

HpCu.X)  =JJ  ds(  - ) exp{A(z+ixcosu+iysinu)  } 


Bo 


dx  exp  (ix( Acosu-k^cosB)  }J  dl  (N^-iN^inB) 
L  C(x) 

x  exp{(kQ+A)z  +  iy(Asinu-kQsing) } 

-/  dx  exp(iAxcosu)/  dl  Aiji  (N  +iN  sinu) 

*L  JC(x)  J  1 

x  exp(Az  +  iAysinu)  . 


NUMERICAL  PROCEDURE 

As  mentioned  earlier,  Frank  (1967)  solved  the  two-dimensional  problem 
where  the  logarithmic  sources  of  Equation  (3.18)  were  distributed  over 
the  hull  cross  section.  Using  the  method  of  linear  superposition,  one 
can  express  the  potential  V  by 

iKP)  -f  o( Q)G(P,Q)dl(Q)  (3.26) 

C 

with  the  unknown  source  strength  a( Q) . 

In  order  to  solve  Equation  (3.26)  two  assumptions  are  introduced 
in  the  numerical  method.  As  shown  in  Figure  8  a  hull  cross  section  is 
described  by  n  offsets  where  dl^  is  the  arclength  between  the  j  and 
j+1  points. 
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Tj-1 

/ 


J<J,rV 


.S'  dli’°3 

J+l<ViV 


Figure  8  -  Cross  Section  of  Hull 


The  first  assumption  is  that  the  source  strength  o(l)  varies  slowly 

enough  to  be  considered  as  constant  0,  over  a  given  arclength  dl .  .  The 

•J  J 

second  assumption  is  that  arc  dl^  can  be  substituted  by  a  straight  line 
dSj  such  that  the  ship  section  could  be  approximated  by  a  chain  of 
straight  lines. 

With  these  assumptions  and  body  boundary  condition  (3.10),  one 
can  obtain  the  following  expression: 


15  -,lM  “MJgCPj.QWUQ) 

J  dl 


V  a  |r/*  G(P.,Q)dl(Q) 
j  =  l  J  mJdl  1 


E  ai  f  G(P  ,Q)ds(Q) 
J-l  j  Jdsi 


(3.27) 
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where  CL  ts  the  constant  source  strength  between  points  j  and  j+1, 

P  the  midpoint  of  the  i-th  arc.  Using  matrix  notation,  we  may  have 


n-1 

f  =  7,  A..O.  , 

1  j.i  lj  J 


i=l,2,  ,n-l 


(3.28) 


where  f  ^  =  — j  and  A  =  ■|jj/dg^G(Pi,Q)ds(Q) .  The  advantage  of 

Equation  (3.27)  is  that  the  term  /G  dl  can  be  evaluated  in  a  closed 
form.  By  increasing  the  number  of  offset  points  this  approximation 
approaches  the  exact  solution.  However  numerical  computations  show 
that  a  relatively  small  number  of  points,  for  example  15  points  for 
a  half  circle,  gives  fairly  good  agreement  with  the  exact  solution. 

We  also  may  increase  the  accuracy  of  solution  either  by  assuming  linearly 
varying  source  strength  over  the  line  segments  ds^  for  a  fixed  number 
of  offset  points,  or  by  integrating  along  an  arc  dl^  instead  of  ds^ 
with  constant  or  linear  ck.  But  the  increased  numerical  complexity 
may  offset  its  merit. 

From  Equations  (3.27)  and  (3.28)  we  can  note  that  for  given  frequency 
co  and  contour  C(x)  the  coefficient  matrix  A_  of  the  forced  motion  or 
diffraction  problem  is  the  same.  Therefore  once  we  obtain  the  inverse 
of  matrix  A^ ,  forced  motion  and  diffraction  can  be  solved  simultaneously. 
This  is  the  consequence  of  replacing  the  Helmoltz  equation  by  the  two- 
dimensional  Laplace's  equation. 


THE  SECOND-ORDER  STEADY  FORCES 

Before  calculating  the  forces,  let  us  discuss  the  numerical  procedure 
for  the  Kochin  function  first.  The  general  form  of  the  Kochin  function  H 
may  be  written  as: 

H(A's)  =j' dx  exp(iA^x  t  dl  f (x;y,z)exp( X2z+i\3y) ,  (3.29) 

L  *C(x) 

where  L  is  the  ship  length  and  f(xjy,z)  has  x  as  a  parameter.  If  we 
assume  that  f  (x;y,z)  varies  smoothly  over  C(x),  then  for  the  contour 
integral  we  might  adopt  the  same  assumptions  as  Frank  did.  That  is, 
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\ 


/dl  f(x;y,z)exp(A  z+iA  y) 
/.  _  _  \  ^  J 


C(x) 


n_1  -  -  r 

X  f(x;y.,z  )/  ds  exp(A  z+iA  y). 
j=l  3  JJdSj 


(3.30) 


Given  Ng  stations  along  ship  length,  we  may  calculate  the  integral  (3.30) 
at  x=Xj  and  denote  its  value  as  F^jA^X^)  and  from  Equation  (3.29) 
we  get, 


H(  A 


i's)  =J  dx  exp(  iA^x)  S1  (x)  ,  for  i=l  ,2 ,  •  •  ’  ,Ng.  (3.31) 


where  S^(X)  is  the  function  to  be  obtained  by  curve  fitting  F  (X^;  A..  ,  Aj) 

Now  it  is  plausible  to  find  a  method  for  obtaining  a  smooth  representation 
for  the  discrete  data  F(X^;A^,A2)  by  the  use  of  the  spline  function  method 
of  curve  fitting.  By  using  this  method  we  may  define  the  interpolating 
spline  as 

S.(x)  =  a.(x-X.)3+  b1(x-X.)2+  c.(x-X.)  +  d.,  (3.32) 


in  the  interval  X^  <.  x  <  X^+^.  The  coefficients  of  the  cubic  polynomial 
are  expressed  in  terms  of  ^ ; A^  ,  A^)  ,  ^ ^i+1’^1 ’^2^  ’  an<3  tlie  seconc* 
derivatives  S^(X±)  ,S'^ (Xi+1)  (see  Appendix  A).  Consequently  we  may  rewrite 
the  Equation(3. 31)  as  follows: 


H(A's) 


Ns"1  i+1 


"s  *  r 

X  /  dx  exp(iA  x)S  (x) 
i=lJ  1  1 


(3.33) 


This  integration  can  be  performed  exactly. 

Equat ion(2.9)  shows  that  the  second-order  steady  force  consists  of 

two  terms,  <AF  >  and  <AF  >  .  The  computation  of  <AF  >  is  straight- 
Id  Bd  IB 

forward  for  a  given  u=b  and  A=k0.  Qn  the  other  hand,  <AFBg>  in  Equation 
(2,11)  is  rather  complicated.  The  reason  is  that  the  Kochin  function  is  a 
function  of  furthermore,  A^A^  t'iemselves  are  function  of  u, 

Figure  2  shows  that  the  A^  goes  to  infinity  as  u  approaches  to  tt/2. 
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The  evaluation  of  the  Kochin  function  becomes  more  difficult  as  A  becomes 
large,  for  the  integrand  of  the  Kochin  function  oscillates  very  rapidly 
along  the  x  and  y  axes  for  large  values  of  A  .  In  the  computation  we 
first  decide  the  angle  u£  which  satisfies  the  following (see  Appendix  B) : 

snnyn 

'  -it/ 2  nr/2  *TT+U  ' 


where  D(u,A  )  = 


u/2-u£  tt-u  3 

■  (/  7  7 

-tt/2+u£  tt/2+u£  TT+ 
A?(u) 

- - - yjy  iHCn+u.Aj 

(1+Axcosu) 


3tt/2-u  f 


du  D(u,A^) 


,  ,  ,  .  1 2  /cosul 

(HCti-Ki,  Ax>  1  (s.nu). 


By  doing  this  we  may  partly  eliminate  the  difficulty  Involved  in 
the  force  integral.  The  contribution  of  each  mode  of  ship  motion  and  the 
diffraction  part  to  the  force  is  examined  separately.  The  numerical 
results  comparing  their  relative  magnitude  (see  Figure  13)  shows  that 
the  major  source  of  the  forces  are  from  heave  and  pitch. 

The  computer  program  based  on  Lin  and  Reed's  theory  has  been  developed 
by  Reed  and  Hubble(1980)  originally.  Extensive  debugging  and  modification 
that  includes  the  irregular  frequency  in  ship  motion  problem,  the  diffrac¬ 
tion  problem  and  the  force  integral  are  made  by  author.  All  the  methods 
has  been  synthesized  into  single  computer  program  by  essentially  combining 
the  ship  motion  program  and  the  Kochin  function  evaluation.  The  ship  motion 
program  consists  of  several  links,  and  the  Kochin  function  is  the  last  link 
to  the  ship  motion  program. 

NUMERICAL  RESULTS 

It  is  obvious  that  computational  accuracy  increases  with  the  number 
of  elements  used  to  approximate  the  body  surface,  and,  in  the  meanwhile, that 
the  cost  for  a  solution  depends  very  strongly  on  the  number  of  points  used. 
The  source  points  should  be  distributed  in  such  a  way  that  the  best  results 
will  be  obtained  with  the  fewest  possible  points.  Naturally,  points 
should  be  concentrated  in  regions  where  the  flow  is  expected  to  change 
rapidly.  In  order  to  decrease  the  computing  cost  within  the  tolerable 


limit  of  accuracy  we  took  the  following  steps;  first,  we  tried  to  bbtain 
the  optimal  minimal  number  of  offset  points  on  the  body.  Second,  within 
the  range  of  frequency  of  practical  interest  appropriate  intervals  of 
frequencies  are  taken  at  which  the  potential,  added-mass  coefficient, 
damping  coefficient,  and  etc.  are  evaluated,  and  the  linear  interpolation 
method  is  used  to  get  those  values  at  the  frequencies  between  the  initially 
chosen  frequencies.  Lastly,  the  integral  of  equation  (2.10)  was  also 
approximated  by  a  finite  sum  of  discretized  integrals. 

The  computations  were  mostly  carried  out  on  the  GDC  6600,  6700 
computer  at  DTNSRDC.  The  evaluation  of  the  influence  coefficient 
matrix  for  the  singularity  strength  and  the  evaluation  of  the  Kochin 
function  were  the  most  time  consuming  parts  of  these  computations. 


Numerical  Examples 


To  facilitate  the  comparison,  the  added  resistance  and  the  lateral 
drift  force  were  non-dimens ionalized  as  follows: 


AF 

x 


added  resistance 
2  2 

PgA  B  /L 


(4.1) 


AF 

y 


lateral  drift  force 
2 

PgA  L 


(4.2) 


Figure  9  shows  the  added-resistance  prediction  for  a  Mariner  hull  form 

at  a  speed  of  F^  =  0.194  with  three  different  headings ,  '$“120° ,  150°  and 

180°  (Note  that  8=  180°  for  head  seas).  A  striking  faet  to  note  is  that 

the  extreme  sensitivity  of  the  added  resistance  to  the  heading  angles. 

For  instance,  at  A /L  =  0.6  which  corresponds  to  approximately  300  feet 

wave  length,  the  added  resistance  can  increase  about  6  times  when  the 

wave  heading  angle  is  changed  from  the  30°  bow  to  60°  bow.  The  results 

in  Figure  9  reveal  that  the  usual  notion  that  the  added  resistance  is 

greater  in  head  seas  is  not  necessarily  true.  In  Figure  10  comparisons  were 

made  between  Salvesen's  calculation  and  <F,  >  at  F  =0.194.  The  major 

11B  n 
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differences  between  Salvesen  (1978)  and  the  present  theory  is  that  Slavesen 
ignored  <Af1bb>  in  Equation(2. 9)  by  assuming  iJ'j.  Relatively  good 

agreement  in  magnitude  is  observed,  but  there  still  exists  discrepancy 
between  them.  There  might  be  two  reasons  for  this;  first,  in  developing 
the  ship  motion  program  some  modifications  were  made  on  the  Frank’s  close- 
fit  method  as  mentioned  earlier.  Secondly,  the  computation  of  the  diffrac¬ 
tion  part  is  quite  different.  In  Figure  11  we  note  that  <AF1t>  compared 
with  <AF^jg>  is  not  small  enough  to  be  neglected.  Figure  12  displays  the 
contribution  of  the  forced  motion  and  the  diffraction  part  to  the  force 

<AF  >  separately  for  the  three  different  headings  at  F  =  0.194. 

1IB  n 

As  would  be  expected,  it  is  seen  that  the  effects  of  the  forced  motion 

decrease  as  X/L  becomes  large,  and  the  maximum  occurs  at  shorter  wavelength 

for  decreasing  heading  angles.  Meanwhile,  the  diffraction  part  seems  to 

act  differently  to  that  of  the  forced  motion. 

Figure  13  shows  relative  magnitudes  of  each  mode  of  ship  motion  for 

<AF  >  separately.  Pitch  and  heave  are  dominant  over  the  others.  It 
1  IB 

is  interesting  to  note  that  the  peak  values  of  pitch  and  heave  for 

8  =  150°  are  slightly  greater  than  those  for  8=  180°.  These  effects  are 

reflected  in  Figure  12  where  the  peak  value  of  the  added  resistance  for 

8  =  150°  is  slightly  greater  than  those  for  8=  180°.  <AF,.>  is  presented 

in  Figures  14  and  15.  Basically  <^F1BB>  consists  of  two  integrals,  i.e., 

X^,  X^  and  the  corresponding  Kochin  functions  H(u,X^),  HCujX^)  respectively. 

As  expected,  Figure  14  shows  that  we  might  neglect  the  highly  oscillatory 

X,  integral  in  the  computation.  In  Figure  15,  <AF.  >  is  expressed  in 

terms  of  each  ship-motion  mode,  diffraction  and  their  interactions. 

The  usefulness  of  any  theory  cannot  be  judged  until  its  prediction 

have  been  compared  with  empirical  data.  Unfortunately,  it  is  difficult 

to  find  experimental  data,  especially  in  an  oblique  seaway.  We  selected 

one  of  the  available  experimental  data;  Series  60,  C  =0.60  at  F  =0.283 

B  n 

in  head  seas.  In  Figure  16  the  theoretical  predictions  of  three  different 


methods  are  presented  together  wi  wo  sets  of  experimental  data  obtained 
by  Sibul  (1971)  and  Strom-Tejsen  et  al.  (1973)  for  this  particular  case. 


In  the  lower  frequency  ranges  not  only  do  all  three  numerical  predictions 
agree  well  with  each  other,  but  also  they  show  fairly  good  agreement 
with  the  experimental  values.  All  the  numerical  values  overpredict  the 
maximum  added  resistances,  however.  Gerritsma  and  Beukelman(1972)  overpre¬ 
dict  it  by  about  a  factor  of  two,  the  present  theory  by  nearly  70%  and 
Salvesen's  method  by  approximately  35%,  In  the  higher  frequency  ranges 
the  present  theory  gives  a  little  better  prediction  when  compared  to  the 
others.  It  is  probably  because  our  computation  of  the  diffraction  poten¬ 
tial,  which  is  the  major  contribution  to  the  added  resistance  in  the  higher 
frequency  ranges,  is  better  than  others.  Two  interesting  things  are  obser¬ 
ved  for  <AF,T>;  first,  <AF,  >  for  this  case  show  good  agreement  with 

1  I  B  1  1  D 

experimental  data.  Secondly  <AF^j.g>  is  slightly  higher  than  the  added- 
resistance  of  Salvesen  for  a  Mariner  but  much  less  than  that  of  Salvesen 
for  Series  60.  We  cannot  give  any  specific  reasons  for  the  descrepancies. 
Only  by  comparing  both  methods  term  by  term,  we  may  find  out  the  differen¬ 
ces. 

The  lateral  drift  forces  for  a  Mariner  are  also  presented  in  Figures 
17  and  18.  Figure  17  shows  the  lateral  drift  force,  with  three 

different  headings  0=  90°,  120°,  150°.  As  would  be  expected,  the  coeffi¬ 
cient  approaches  0,5  as  the  frequency  becomes  higher  for  the  beam-sea 
case,  but  the  predictions  give  small  negative  values  for  oblique  seas  in 

the  higher  frequency  ranges,  <AF  >  and  '  AF  >  are  plotted  separately 

Z  1  D  4UD 

in  Figure  18.  The  relative  magnitude  of  <AF  >  compared  with  <AF„T  > 

Zc> d  Z  Id 

is  much  smaller  than  that  of  <AF.  >  with  <AF,  >.  Comparisons  of  added 

IBB  1  IB 

resistance  and  lateral  drift  force  with  experimental  values  will  be  given 
in  the  future. 

CONCLUDING  REMARKS 

The  second-order  steady  forces  have  been  considered  for  a  ship  in 
regular  waves  of  arbitrary  headings.  The  accurate  prediction  of  these 
forces  is  of  considerable  importance  for  estimating  the  powering  require¬ 
ment  in  waves,  assessment  of  seakeeping  qualities,  and  the  position 
control  of  ships.  Based  on  Lin  and  Reed  (1976)  a  new  numerical  scheme 
has  been  developed  for  predicting  the  added  resistance  and  drift  force. 
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The  summary  of  the  findings  in  this  study  is  as  follows: 

(1)  It  is  not  so  easy  to  verify  the  computed  second-order  steady 
forces  satisfactorily  because  of  the  lack  of  reliable  experimental  data. 

(2)  The  effect  of  wave  heading  on  added  resistance  is  significant. 
The  maximum  added  resistance  does  not  necessarily  occur  in  the  head  seas 
only. 

(3)  For  decreasing  heading  angle  or  X/L,  the  contribution  of  the 
motion  to  the  added  resistance  decreases  while  that  of  the  diffraction 
part  increases; 

(4)  There  exists  a  contradiction  between  the  assumptions  made  in 

ship  motion  and  the  diffraction  problem.  In  the  ship  motion  problem  the 

justification  of  the  strip  theory  is  made  by  assuming  A=0(B) ,  i.e.,  the 

wavelength  is  approximately  of  the  order  of  B,  while  in  diffraction 
;  2  2 

problems  the  term  (k  cos  8)  is  dropped  in  order  to  replace  the  Helmholtz 
equation  by  Laplace's  equation.  That  means  that  in  the  diffraction 
problem  the  assumption  of  X^B  is  made. 

It  is  fair  to  say  that  the  theory  of  Lin  and  Reed  is  mathematically 
sound,  but  there  still  exists  a  gap  between  the  theory  and  its  practical 
applicability.  In  concluding  this  work  we  like  to  make  some  suggestions 
for  future  study.  First  of  all,  to  develop  a  numerical  prediction 
method  of  the  second-order  steady  forces  it  is  desired  to  use  the 
most  accurate  method  for  predicting  ship  motion.  For  the  diffraction 
problem  we  have  to  solve  either  three-dimensional  Laplace's  equation 
or  Helmholtz  equation  in  high  frequency  ranges. 
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APPENDIX  A 


Spline  Curve  Fitting 

Suppose  we  wish  to  approximate  a  continuous  and  differentiable  function 
F(x)  on  the  interval  {0,L}  in  a  piecewise  fashion,  using  low-degree 
interpolating  polynomials  over  nonoverlapping  subintervals  of  {0,L}. 

Let  the  base  points  be  0  =  x^<  x^<...  xn  =  L,  the  corresponding 

functional  vaules  be  y^=F(x^),  i=0, 1 , 2, . . . ,n,  and  interpolating  function 
for  {0,L}  be  S^(x).  We  shall  require  that  S^(x)  be  continuous  on  {0,L} 
and  possess  continuous  first  and  second  derivatives  for  all  x  in  {0,L>. 


Let  S^(x)  coincide  with  a  third-degree  polynomial  on  each  interval , i. e. , 
S.(x)  -  »1(x-X1)3  +  b^x-X.)2  +  c1(x-X1)  +  dj,  (A.  1 ) 


<  Xt<  x<  X1+1  ) 


then 


S^Cx)  =  Sa^Cx-X^)  +  2b^(x-X  J  +  c ^ 

(A.  2) 

and 

S^x)  =  6a^(x-X^)  +  2bi> 

(A.  3) 

Now,  by  setting 

Si(Xt)  =  Yi  and  S  (Xi+1)  «  Y.+1>  Equation(A. 1) 

yields 

Y.  =  d., 

1  X* 

(A.  4) 

Yi+l  =  aih3i+bih2i  +  cihi  +  di’ 

(A- 5) 

where  h  =X. , ,  - 
i  l+I 

xr 

At  each  of  the 

interior  points,  we  set 

s;««>  ■  si+lixi+1), 

or 

3aihi  +  2b!hi  +  ci  '  ci+l. 

(A.  6) 

and 

Si<XW>  *  S’Ll(Xl+l)- 

or 

6a, h.  +  2b.  =  2b  ,  for  i=l ,2, . . . ,n-2. 

i  i  l  l+l 

(A.  7) 

Since  the  second  derivative  is  a  piecewise  linear  function  of  x, 
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where  X^ 
Equat ions 


S'^(x)  -  X  +  - -i-  (x_x  ) 

i+1  "  Xi  1 


=  S’^X^  and  Xi+1 


Si(Xi+l>  "  Si+l(Xi+l>- 


(A. 3)  and  (A. 8)  there  results 


2bi  =  Xi’ 

6ai  “  (  Ai+rAi)/hi' 


(A.  8) 


Thus,  comparing 


(A.  9) 
(A.  10) 


From  Equations  (A. 4) , (A. 5) , (A. 9)  and  (A. 10)  the  coefficients  of  the  cubic 
polynomial  in  the  interval  {X^.X^}  can  be  expressed  in  terms  of  Y  ,  Y  + 
and  X^,  X^+^  as  follows: 


ai =  (Ai+rAi)/6hi‘ 


bi  -  V2- 


ci =  (Yi+rYi)/hi  '  (2Ai+Ai+i  )/6* 


(A.  11) 


dt  =  Yr 


Substituting  Equat ion(A. 11)  into  Equation(A. 6)  and  after  some  algebric 
manipulation,  we  obtain  the  basic  equation  of  Spline  technique  as  follows: 


h. 

l 

hi+l 

_2<hi+hi+i» 

Ai  + 

Ai+1  + 

2<hi+Vi> 

3 

Y  -Y 

Y 

_  _ 

i+2 

i+l  *i+l 

i+2 


hi+hi+l 


i+1 


(A.  12) 


for  i=l ,2, . . . ,n-l ,n-2. 
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APPENDIX  B 


Let  us  denote 


X  Integral  of  <AF  > 

1  DD 


G(X1,u) 


X^(u)  cos  u 
(1  +  4tcos  u) 


%  |  H(-n+u,X1) 


(B.l) 


71/2  TT-U0  3TT/2 

I  ={  /  -J  -J  }G (Xj  ,u)du , 

—IT/ 2  n/2  "-Hi 


where  T=u.'U/g,  H(tt+u,X^)  and  X^(u)  are  defined  in  Equations  (3,1)  and  (3.2), 
respectively.  As  shown  in  Figure  2,  X^  becomes  infinite  as  u  approaches 
to  tt/2,  and  if  becomes  large,  the  integrand  of  Equat  ion (3 . 1 )  oscillates 
so  rapidly  that  it  is  difficult  to  evaluate  the  integral  properly. 

We  will  examine  the  second  integral  of  F.quation(B,  2)  first.  Let  us  assume 
that  there  exist  a  small  positive  angle  u^  such  that  as  u^  0,  we  have 


O  uo 

J"  C(X,u)du  -  J"  G(X,u)du, 


n/2  ti/2-Ki 

where  the  subscript  "1"  of  X  is  omitted  for  the  brevity  sake.  Let  us 
define  the  difference  between  the  exact  and  the  approximate  value  by 


Error 


r  tt/2+u 

2/  £ 

‘'tt/2 


G  (  X ,y)  du, 


(B.4) 


and  determine  the  angle  such  that  the  Error  is  within  an  acceptable 
limit.  In  order  to  do  this,  let  us  examine  the  magnitude  of  the  Kochin 
f unct iom 


H(tt+u,X) 


*  f  dx  exp(iXxcos  u) f  dlX^(n.+in  sin  u) 

-1  */C(x) 

*  exp{X(z  +  iysin  u)}  ,  for  v/2£y  <tt/2+u  .  (B.5) 


Here  we  neglect  the  term  (iun)  in  Equation(3.1)  under  the  assumption  that 
X  is  very  large.  Using  Equation(3. 30)  we  can  approximate  the  Kochin 
function  for  a  large  value  of  X  as  follows? 

-1  n-1 

H(tt-Ki,X)  =  /  dx  exp(iXxcos  u)  Z  f  (x;,y.,z  ,u) 

J-1  j-1  J  J  J 

x  J"  ds  Xexp{ X (z+iysin  u)} 


where  f  (x;  y .  ,  z .  ,u)=  (n  -tin.  sin  u)  and  (y  ,z.)  is  the  magnitude  of 
J  3  J  J  -^J  J  J 

the  ith  segment,  and  furthermore,  considering  the  fact  that  | exp ( iAxcosu)  |  , 

|  exp(lXysinu)  |<.  1.0,  we  may  have 


/  n-1  _  ij- 

I H  (iT+u ,  X )  j  <  I  /  dx  Z  f  (x;y,,z  ,u )/ 

-1  j=l  J  J  J  z. 


dz  Xexp (Xz) 


1  exp(-X|z|  )  |D(u) 


(B.6) 


1  n-1 


/X  u  —  X 

dx  j  =1  f  j  (x;yj  >u) 


(B.7) 


and  z  .  is  the  minimum  value  of  z.  ,-z,  of  all  the  cross  sections. 
'  min  '  j+1  j  1 

Combining  Equations  (B.4)  and  (B.6),  and  using  the  mean-value  theorem, 

we  have 


tt/2-Hi 


<  f  C  du  -V.(u).cos_u  r..  exp{-2>|z|  )|n(u)|? 

(1  +  4  TC  o  su ) J  min 


V 


\ 


|D(u)|  2f  '  du 

J  n t  /  O 


A  (u)cOS  U 

(1  +  4tcosu) 


5“  exp{-2A|z|min},  (B.8) 


where  D(u)  defined  in  Equation(B, 7)  is  a  smooth  function  of  u  and  u  is  the 
value  between  tt/2  and  tt/2+u£.  By  change  of  variable  the  integral  in 
Equation(B,8)  will  become 


'  ue  A  (u) (-sin  u) 

du  — — — - - — 


~ -  ~  exp{-2X  (z  lmin>- 

(1  -  4xsmu) 


(B ,  9) 


Assuming  u£  is  small  and  keeping  the  leading  term  only,  we  obtain  the 
following! 


/  "  eXp("TI|z,niin)  . 

J  4  3  t  u 

0  Xu 


(B. 10) 


We  will  denote  Equat ion (B. 10)  as  A  and  again  apply  the  change  of  variable 
u2  =  v,  the  result  will  be 
2 

A  =  |  -ex-P..^/ALl  dv  =  a/(28)  exp(-6/u2),  (B.ll) 

“  J  0  v 

where  i  =  v2/x  and  8  =  2v|z|  /  2, 

min  x 

Let  us  put 


1/(28)  exp(-P/u2)  10"P, 


(B.12) 


where  the  arbitrary  positive  value  P  will  decide  the  accuracy  of  compu¬ 
tation  and  the  angle  u  simultaneously,  Equat ion(B. 12)  can  be  rewritten 


~  10P  <_exp(6/u£), 

or 

ln(  fg  10P  )  <  6/u£  , 

and  if  yg  10P  _>  1.0  ,  then  we  will  have 

u£  <{3/ln(al0P/ (26))}  '"2  .  (B.13) 

This  value  u£  gives  the  absolute  magnitude  of  the  Integral  of  Equation 
(B.4),  not  a  relative  size  of  error,  i.e,,  the  ratio  of  the  magnitude  of 
Equation(B, 4)  to  that  of  Equation(B, 3) .  Because  of  the  difficulty  of 
integration  of  Equat ion(B. 3) ,  we  cannot  confirm  the  validity  of  Equation 
(B.13)  but,  in  the  following  example,  we  show  the  numerical  vlaues  of 
Equation(B, 3)  for  several  different  P’s.  In  stead  of  a  actual  ship,  a 
rectangular  barge  is  considered  for  an  example  with  the  following  condi¬ 
tions* 


L/B  =  7.5, 

B/T  =  2.0, 

F  =0.2, 
n 


6  =  135  (degree) , 

A/L  =  2.0, 

( z  |  .  /B  =  0.1. 

1  'min 


In  addition  to  these  the  velocity  potential  in  Equation(B.  5)  is  assumed 
to  be  constant.  The  Kochin  function  H(n+u,\)  is  obtained  in  the  close 
form, i. e. , 


H(tt-Ki,A)  = 


4sln(Alsln  u) ,  e 


-A 


A  cos  u 


{ 


sAr\C\sAn,  lO  _  isinucos ( Acos  u)(e  ^-1)  } 
sin  u 


where  l=L/2,  and  the  integral  of  Equat ion (B. 3)  may  be  summarized  as  follows; 
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p 

u 

e 

(degree) 

A(u  ) 

€ 

/* ''  o 

I  = f  G  (X ,u)du 

tt/2+u^ 

Vh 

3 

25.68 

4.87 

-9043.66 

1.0000 

4 

22.51 

6.98 

-9073.47 

1.0033 

5 

20.28 

9.11 

-9114.01 

1.0078 

6 

18.61 

11.26 

-9144.22 

1.0111 

7 

17.29 

13.43 

-9169.18 

1.0139 

8 

16.22 

16.62 

-9187.22 

1.0159 

9 

15.32 

17.83 

-9211.14 

1 . 0185 

10 

14.56 

20.05 

-9225.95 

1 . 02  02 

where  u0=cos  ^ (X/ At) =7 1 . 39° .  Though  P  increases  f i on  3  to  10,  the  total 
increment  of  Equat ion (B. 3)  for  P=10  is  about  2%  of  the  integration  for  P=3. 
The  integral  does  not  converge  as  fast  as  P  grows,  but  it  does  approach 
to  finite  value  as  P  becomes  large. 

The  third  integral  of  Equat ion (B. 2 )  is  identical  to  the  second  one,  and 
because  X  ( rr  /  2  — u  ^ )  >  X  (it/ 2+u.  )  ,  the  angle  u£  determined  in  Equat  ion  (B .  1  3) 
can  also  be  used  for  the  first  integral  of  Equat ion (B, 2) .  Consequently, 
we  have 


I 


-it/  2+u 

c 


"/  2+u 


G(‘ ,u)du. 


(B . 14 ) 


SI 
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Figure  5  -  Added-Mass,  Damping  Coefficients  and  Heave,  Pitch  Amplitude 
for  Mariner  at  F  =0.20,  ('= '  1 5 n 
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Figure  14  -  Comparison  between  X  and  X  Integral  of  <AF_ > 
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at  F  »0.194 
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Added  Resistance/ (pgA  B  / L) 


Figure  16  -  Comparison  of  Experimental  and  Theoretical  Added  Resistance 
for  Series  60  Hull  with  C  =0. 60  at  F  -0.283,  6-180° 
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Figure  17  -  Lateral  Drift  Force,  <AF->,  for  Mariner  at  F  -0.194 
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DTNSRDC  ISSUES  THREE  TYPES  OF  REPORTS 

1.  DTNSRDC  REPORTS.  A  FORMAL  SERIES,  CONTAIN  INFORMATION  OF  PERMANENT  *!CH 
NICAL  VALUE.  THEY  CARRY  A  CONSECUTIVE  NUMERICAL  IDENTIFICATION  REGARDL  OF 
THEIR  CLASSIFICATION  OR  THE  ORIGINATING  DEPARTMENT. 

2.  DEPARTMENTAL  REPORTS,  A  SEMIFORMAL  SERIES,  CONTAIN  INFORMATIOl  OF  A  PRLlIW 
INARY,  TEMPORARY,  OR  PROPRIETARY  NATURE  OR  OF  LIMITED  INTEREST  OR  SIGNIFICANCE. 
THEY  CARRY  A  DEPARTMENTAL  ALPHANUMERICAL  IDENTIFICATION. 

3.  TECHNICAL  MEMORANDA,  AN  INFORMAL  SERIES,  CONTAIN  TECHNICAL  DOCUMENTATION 
OF  LIMITED  USE  AND  IN  I  EREST.  THEY  ARE  PRIMARILY  WORKING  PAPERS  INTENDED  FOR  IN 
TERNAL  USE.  THEY  CARRY  AN  IDENTIFYING  NUMBER  WHICH  INDICATES  THEIR  TYPE  AND  THE 
NUMFRICAL  CODE  OF  THF  ORIGINATING  DEPARTMENT.  ANY  DISTRIBUTION  OUTSIDE  DTNSRDC 
MUST  BE  APPROVED  BY  THE  HEAD  OF  THE  ORIGINATING  DEPARTMENT  ON  A  CASE  BY  CASE 
BASIS. 
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